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Abstract 

The  performance  of  polymer  electrolyte  membrane  (PEM)  fuel  cell  stacks  can  be  improved  significantly  by  optimizing  the  design  and 
operating  conditions.  In  this  study,  performance  modeling  and  optimization  of  a  PEM  fuel  cell  stack  have  been  conducted.  The  pressure  and 
molar  flow  rate  distributions  for  the  fuel  and  oxidant  streams  in  the  stack  are  determined  with  a  flow  network  model  incorporating  the  minor 
losses.  The  distributions  are  then  used  in  the  single  cell  model  developed  previously  to  evaluate  the  performance  of  PEM  fuel  cell  stack. 
Analysis  has  been  carried  out  for  different  fuel  and  oxidant  flow  configurations  and  bipolar  plate  designs.  It  was  found  that  the  minor  losses 
increase  the  stack  operating  pressure  and  the  power  requirement  for  oxidant  supply  and  alter  the  cell-to-cell  voltage  variations  in  the  stack.  A 
symmetric  double  inlet-single  outlet  topology  provides  optimal  stack  performance  with  reasonably  low  compressor  power  requirement  for 
the  reactant  flow  and  minimum  cell-to-cell  voltage  variations.  The  stack  performance  is  considerably  affected  by  the  size  and  the  number  of 
flow  channels  on  bipolar  plate.  Optimal  stack  performance  requires  the  matching  of  the  stack  manifold  designs,  flow  channels  on  the  bipolar 
plates  and  the  stack  operating  conditions. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cells  convert 
the  chemical  energy  of  hydrogen  and  oxygen  directly  and 
efficiently  into  electrical  energy  and  are  widely  regarded  as 
alternative  stationary  and  mobile  power  source.  The  main 
characteristics  of  PEM  fuel  cells  are:  they  produce  water  as  a 
byproduct;  they  have  higher  efficiency  when  compared  with 
heat  engines;  they  operate  at  low  temperatures  (up  to  90  °C), 
which  allows  a  fast  start-up;  they  use  a  solid  polymer  as  the 
electrolyte,  which  reduces  concerns  related  to  construction, 
transportation,  and  safety.  Before  PEM  fuel  cells  can  be 
successfully  commercialized,  the  production  cost  must 
be  reduced  from  the  current  estimate  of  approximately  $ 
200/kW  to  $  30/kW  [1]. 
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The  output  voltage  of  a  single  PEM  fuel  cell  is  limited, 
usually  less  than  1  V.  In  order  to  produce  a  useful  voltage 
for  practical  applications,  several  unit  cells  are  connected  in 
series  to  form  a  fuel  cell  stack.  The  output  voltage  depends 
on  the  number  of  unit  cells.  Fuel  and  oxidant  flow  through 
the  fuel  cells  in  multiple  channels  arranged  in  a  complex 
flow  network.  This  produces  major  challenges  for  fuel  cell 
designers  since  the  flow  distribution  within  a  fuel  cell  stack 
could  have  a  significant  impact  on  fuel  cell  performance  and 
efficiency,  and  the  heat  and  water  management  strategies, 
which  are  successful  for  single  PEM  fuel  cells,  are  difficult  to 
implement  in  a  stack  environment.  As  a  result,  the  efficiency 
and  power  output  of  a  PEM  fuel  cell  operating  within  a  stack 
are  lower  than  the  performance  of  a  PEM  fuel  cell  operating 
independently  [2]. 

Over  the  past  few  years,  a  significant  amount  of  research 
has  been  devoted  to  the  study  of  PEM  fuel  cell  stacks  and  sev¬ 
eral  empirical  and  mathematical  models  have  been  developed 
for  the  purpose  of  understanding  and  predicting  PEM  fuel 
cell  performance.  A  comprehensive  review  of  the  published 
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Nomenclature 

C  molar  concentration  (mol  m-3) 

Cf  wall  friction  coefficient 

dh  flow  channel  hydraulic  diameter  (m) 

D  diffusion  coefficient  (m2  s-1) 

Dh  manifold  hydraulic  diameter  (m) 

E  cell  or  stack  voltage  (V) 

/pt  mass  ratio  of  platinum  to  carbon  support 
F  Faraday  constant  (96485  C  mol- 1 ) 

H  bipolar  plate  effective  height  (m) 

J  cell  current  density  (A  m-2) 

Ke  electrokinetic  permeability  of  polymer  elec¬ 

trolyte  membrane  (m2) 

Kv  hydraulic  permeability  of  polymer  electrolyte 

membrane  (m2) 

/  flow  channel  length  (m) 

im  fraction  of  catalyst  layer  void  space  occupied 
by  liquid  water 

L  control  volume  length  (m) 

mpt  platinum  mass  loading  per  unit  electrode  area 

(kgm-2) 

m  mass  flow  rate  (kg  s-  ) 

n  exponent  used  to  represent  frictional/minor 

losses 

N  number  of  cells/loops/branches/flow  chan¬ 

nels/turns 

N  molar  flow  rate  (mol  s  ) 

No  initially  assumed  molar  flow  rate  (mol  s-  ) 

Nr  rate  of  reactant  consumption  in  the  catalyst 

layer  (mol  s-1) 

P  pressure  (Pa) 

Pc  compressor  power  (W) 

r  flow  resistance  coefficient 

R  universal  gas  constant  (8.314kJkmol-1  K-1) 

Re  Reynolds  number 

S  stoichiometry 

Se  voltage  spread  (%) 

T  stack  temperature  (K) 

V  average  velocity  (ms-1) 

W  bipolar  plate  effective  width  (m) 

Greek  letters 

a  ratio  of  pressure  loss  in  the  bipolar  flow  chan¬ 

nels  to  that  of  the  manifold 
8  thickness  (m) 

A  difference 

0  overpotential  (V) 

0  flow  direction  convention  (+1  for  clockwise, 

—  1  for  counter  clockwise) 

ks  conductivity  of  catalyst  layer  solid  phase 

(S  m-1) 

fi  viscosity  (N  sm-2) 

p  density  (kgm-3) 


..bulk 

resistivity  of  electrode  backing  (Q  m) 

0 

porosity 

Subscripts 

a 

anode 

bend 

flow  channel  bends 

bP 

bipolar  plate 

branch 

branch  (control  volume) 

c 

cathode/flow  channel/catalyst  layer 

cell 

fuel  cell 

cp 

cooling  plate 

drag 

dragged  molecules 

e 

electrode  backing 

f 

friction 

H+ 

protons 

h2 

hydrogen 

o2 

oxygen 

i 

loop  number 

in,  inlet 

in/inlet  value 

j 

branch  number 

loop 

loop 

max 

maximum  value 

min 

minimum  value 

m 

minor  loss/polymer  electrolyte  membrane 

out,  outlet  out/outlet  value 

r 

reacting 

rev 

reversible 

stack 

stack 

t 

total/turn 

models  can  be  found  in  Baschuk  and  Li  [3,4].  The  majority 
of  these  stack  models  are  extensions  of  the  empirical  and/or 
mathematical  PEM  fuel  cell  models,  originally  developed  for 
a  single  PEM  fuel  cell.  Reactant  distributions  within  the  stack 
is  modeled  by  treating  the  stack  manifold  and  gas  flow  chan¬ 
nels  as  hydraulic  network,  and  the  voltages  of  the  cells  in 
the  stack  are  determined  using  single  cell  models.  Generally, 
these  models  suffer  from  two  limitations:  minor  losses  are  ne¬ 
glected  and  the  stack  performance  is  studied  under  a  specified 
manifold  configuration.  As  will  be  shown  in  this  study,  the 
minor  losses  have  significant  influence  for  the  typical  PEM 
fuel  cell  stack  designs.  As  a  result,  these  model  predictions 
may  not  be  adequately  justified. 

To  this  end,  the  present  study  formulates  a  PEM  fuel  cell 
stack  by  incorporating  multiple  inlet/outlet  topologies  and 
minor  loss  effects.  The  fuel  and  oxidant  flow  distributions 
within  the  stack  are  determined  by  applying  the  conserva¬ 
tion  of  mass  and  energy  equations  to  the  stack  flow  network. 
Cell  voltages  are  determined  using  the  steady  state,  isother¬ 
mal,  single  cell  model  developed  previously  by  Baschuk  and 
Li  [3].  A  variety  of  flow  configurations  for  a  fuel  cell  stack 
operating  with  humidified  hydrogen  and  air  are  simulated, 
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Fig.  1.  Schematic  illustration  of  a  PEM  fuel  cell  stack. 


strategies  for  reducing  the  unequal  distribution  of  reactants 
within  the  stack  are  examined,  methods  for  the  improvement 
of  stack  performance  are  described,  optimization  of  stack  de¬ 
sign  is  performed,  and  optimal  stack  design  is  obtained. 


2.  Model  formulation 


in  order  to  remove  the  heat  produced  by  the  PEM  fuel  cells 
and  maintain  a  constant  stack  temperature. 

The  stack  performance,  often  measured  in  terms  of  the 
stack  voltage,  EstSLCk,  can  be  determined  by: 


A^cell  Ncell 

^stack  —  ^  ^  ^cell  —  ^  ^  Ocp 
1  1 


A  diagram  of  a  typical  section  of  a  fuel  cell  stack  is 
shown  in  Fig.  1.  The  individual  cells,  referred  to  as  mem¬ 
brane  electrode  assemblies  (MEA),  are  composed  of  a  mem¬ 
brane  electrolyte  sandwiched  between  two  porous  electrodes. 
The  MEAs  produce  direct  current  electricity.  Bipolar  plates, 
which  are  electrically  conducting,  separate  the  MEAs  as  well 
as  provide  a  means  for  delivery  of  the  fuel  and  oxidant  to 
the  catalyst  layer  at  the  electrode-membrane  interface.  The 
basic  unit  is  repeated  to  build  up  a  stack.  A  complete  multi¬ 
cell  stack  may  include  cooling  plates,  which  are  specially 
designed,  and  similar  to  bipolar  plates.  These  plates  are  de¬ 
signed  to  provide  a  uniform  stack  temperature. 

The  fuel  and  oxidant  for  each  PEM  fuel  cell  are  supplied 
by  the  stack  manifold,  with  the  anode  manifold  supplying 
fuel  and  the  cathode  manifold  supplying  oxidant.  The  major 
component  of  the  fuel  for  a  PEM  fuel  cell  is  hydrogen,  with 
carbon  dioxide  and  carbon  monoxide  being  present  if  refor¬ 
mate  fuel  is  utilized.  The  oxidant  used  in  a  PEM  fuel  cell 
is  oxygen,  with  nitrogen  being  present  if  air  is  used  as  the 
oxygen  supply.  The  gas  flow  channels  remove  the  water  pro¬ 
duced  by  the  electrochemical  reactions  within  the  membrane 
electrode  assembly  and  supply  the  humidity  required  to  avoid 
polymer  electrolyte  membrane  from  dehydration.  In  addition 
to  fuel  and  oxidant,  water  is  circulated  through  cooling  plates 


where  Aceii  is  the  total  number  of  fuel  cells  in  the  stack,  ECQ\\ 
is  the  voltage  of  each  cell,  and  ijc p  is  the  ohmic  loss  due  to  a 
cooling  plate,  and  is  determined  by  treating  the  ohmic  loss  for 
the  cooling  plate  similar  to  the  ohmic  loss  for  bipolar  plate. 

The  stack  model  presented  here  consists  of  two  parts:  the 
single  cell  model  and  the  stack  flow  network  model.  The  sin¬ 
gle  cell  model  determines  the  voltage  of  each  cell  in  the  stack 
based  on  the  cell  inlet  pressure,  temperature,  stoichiometry, 
and  reactant  composition  in  the  gas  flow  channels,  as  well 
as  the  current  density  and  design  parameters.  The  flow  con¬ 
ditions  around  each  cell  is  determined  using  the  stack  flow 
model. 

2.7.  Single  cell  model 

The  voltage  of  each  cell  in  the  stack  is  determined  using 
the  single  cell  model  developed  by  Baschuk  and  Li  (2002) 
[3].  According  to  this  model,  the  output  voltage  of  a  single 
cell  can  be  determined  as  follows: 

^cell  =  T^rev  ~  Osl  ~  \0c\  2?]e  T]m  (2) 

where  EYev  is  the  reversible  cell  voltage,  r/a  and  r]c  are  the 
overpotentials  attributed  to  the  anode  and  cathode  catalyst 
layers,  respectively.  The  voltage  losses  caused  by  the  bipolar 
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plate,  electrode  backing  and  polymer  electrolyte  membrane 
are  denoted  by  ribp,  and  Om,  respectively. 

The  reversible  cell  voltage  is  the  cell  potential  obtained 
at  thermodynamic  equilibrium.  It  is  a  function  of  tempera¬ 
ture  and  reactant  concentration  through  a  modified  version 
of  the  Nernst  equation.  The  cell  voltage  is  reduced  from  the 
reversible  cell  voltage  by  the  overpotentials  associated  with 
the  various  components  of  the  PEM  fuel  cell.  The  anode  and 
cathode  catalyst  layer  overpotentials  are  found  by  consid¬ 
ering  species  conservation,  and  proton  and  electron  migra¬ 
tion  within  the  catalyst  layers.  Proton  and  electron  migration 
through  the  catalyst  layers  are  related  to  the  protonic  and  elec¬ 
trical  current  using  Ohm’s  law.  Species  conservation  requires 
modeling  of  reaction  kinetics  and  mass  transport.  Oxygen  re¬ 
duction  is  modeled  with  the  Butler- Volmer  equation  in  the 
cathode  catalyst  layer,  while  in  the  anode  catalyst  layer  the 
adsorption  and  desorption  of  H2,  CO  and  O2,  the  electro¬ 
oxidation  of  the  adsorbed  hydrogen  and  carbon  monoxide, 
and  the  heterogeneous  oxidation  of  H2  and  CO  by  O2  are  in¬ 
cluded  in  the  reaction  kinetics.  The  voltage  losses  attributed 
to  the  bipolar  plate  and  electrode  backing  are  the  result  of 
electron  migration;  the  overpotential  is  calculated  by  con¬ 
sidering  the  electrode  backing  and  bipolar  plate  as  electrical 
resistances.  The  overpotential  associated  with  the  proton  mi¬ 
gration  in  the  polymer  electrolyte  membrane  is  determined 
by  the  Nernst-Planck  equation  assuming  a  constant  conduc¬ 
tivity  for  the  fully  hydrated  polymer  electrolyte  membrane. 

Determination  of  the  reversible  cell  voltage  and  overpo¬ 
tentials  requires  several  operating  and  design  parameters.  De¬ 
sign  parameters  depend  on  the  manufacture  of  the  PEM  fuel 
cell  and  include  properties,  such  as  conductivity  and  porosity, 
and  geometric  dimensions.  The  operating  parameters  include 
current  density,  temperature,  pressure,  reactant  composition 
and  stoichiometry. 

Due  to  the  series  connection,  the  current  density  in  each 
cell  will  be  equal.  The  circulation  of  cooling  water  maintains 
a  uniform  temperature  for  each  cell  in  the  stack.  However, 
the  pressure,  reactant  composition  and  stoichiometry  can  vary 
from  cell  to  cell  if  the  mass  flow  rate  and  pressure  distributions 
within  the  stack  are  unequal;  this  variation  is  analyzed  with 
the  stack  flow  model  described  in  the  next  section. 


2.2.  Stack  flow  model 

The  cross-section  of  a  PEM  fuel  cell  stack  was  shown  be¬ 
fore  in  Fig.  1.  Fig.  2  illustrates  the  structure  of  a  typical  bipo¬ 
lar  plate  with  manifolds  and  three  flow  channels  arranged  in 
serpentine  configuration.  To  model  the  flow  distribution  in 
the  stack,  the  complex  fuel  and/or  oxidant  flow  paths  con¬ 
sisting  of  the  stack’s  inlet(s),  outlet,  manifolds  and  gas  flow 
channels  can  be  reduced  into  a  graphical  flow  network  as  de¬ 
picted  in  Fig.  3.  In  this  diagram,  fuel  cells  are  surrounded  by 
the  manifolds  and  flow  channels.  The  top  manifold  supplies 
the  reactants  to  the  gas  flow  channels  of  the  fuel  cell  stack 
and  the  gas  flow  channels  exit  into  the  bottom  manifold. 


Nc  =  3  Nt  =  4 


Fig.  2.  Illustration  of  the  bipolar  plate  with  serpentine  flow  channel  config¬ 
uration  with  three  channels  per  plate  and  four  U-turns  per  channel. 

In  practice,  the  electrochemical  reactions  start  at  the  be¬ 
ginning  of  flow  channels  immediately  after  the  reactants  are 
brought  into  contact  with  the  catalyst  layer.  The  reaction  rates 
vary  along  the  flow  channels  as  the  reactant  composition  and 
pressure  change.  However,  the  present  network  model  uses 
a  control  volume  approach  to  simplify  the  detailed  analysis. 
The  top  and  bottom  manifolds  are  segmented  to  account  for 
the  flow  and  composition  change  as  they  interact  with  the  flow 
channels.  Each  gas  flow  channel  is  divided  into  two  sections 
or  branches  of  identical  lengths.  The  upstream  channel  starts 
from  the  top  manifold  and  extends  half  way  through  the  chan¬ 
nels  length.  The  downstream  channel  starts  from  the  channel 
midpoint  to  the  bottom  manifold.  The  boundary  between  the 
two  sections  is  indicated  by  the  plane  A- A  in  Fig.  2.  It  is 
assumed  that  the  pressure  and  compositions  are  uniform  in 
each  control  volume.  Mass  transfers  due  to  electrochemical 
reactions  are  considered  to  take  place  at  average  operating 
conditions  determined  at  the  middle  of  the  flow  channels. 
These  reaction  sites  are  denoted  by  the  symbol  “(g)”  in  Fig.  3. 
Therefore,  the  simplified  stack  flow  network  consists  of  sev¬ 
eral  loops,  with  each  loop  comprising  six  branches:  two  in 
the  manifolds,  and  two  in  each  flow  channels.  The  interfaces 
at  which  manifold  branches  meet  are  represented  by  nodes 
as  illustrated  in  Fig.  3. 

With  the  present  multiple  inlet  stack  model,  numerous 
flow  configurations  can  be  examined  for  optimum  PEM  fuel 
cell  stack  performance.  The  inlets  supply  humidified  gaseous 
hydrogen  and  air  (oxygen)  to  the  anode  and  cathode  sides, 
while  the  outlets  remove  the  reaction  product  and  excess  re¬ 
actants. 

2.2.7.  Conservation  of  mass 

The  inlet  hydrogen  and  oxygen  molar  flow  rates  to  the 
anode  and  cathode  sides  of  the  stack  are  determined  by 
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Fig.  3.  Simplified  flow  network  configuration  of  a  PEM  fuel  cell. 


Attack  _  ^a^cell/AceH  Xrstack  _  ^c^cell^cell 

nh2  -  ^  ^ 

where  <Sa  and  Sc  are  the  stack  anode  and  cathode  stoichiome¬ 
tries,  respectively.  Total  inlet  molar  flow  rates  can  be  calcu¬ 
lated  by  adding  the  water  vapor  and  nitrogen  to  the  inlet 
streams.  The  maximum  amount  of  water  vapor  coming  into 
the  stack  with  the  reactants  corresponds  to  100%  relative  hu¬ 
midity  in  the  incoming  gas  streams. 

The  distribution  of  incoming  components  is  governed  by 
the  conservation  laws.  The  conservation  of  mass  principle 
requires  that  total  mass  flow  into  and  out  of  each  node  in  the 
flow  network  (Fig.  3)  must  be  equal,  or 

mt  =  mt  (for  all  nodes)  (3) 

in  out 


The  exit  stream  from  the  anode  is  assumed  to  be  fully  sat¬ 
urated  with  no  condensed  water.  Since  H2  is  consumed  in 
the  anode  reaction  sites,  the  excess  water  molecules  are  con¬ 
sidered  to  migrate  through  the  polymer  electrolyte  membrane 
towards  the  cathode  side  as  a  net  result  of  electroosmotic  drag 
and  back  diffusion.  The  molar  consumption  in  the  anode  can 
be  expressed  as: 

A  TV  Anode  =  A^H2  +  ^H20  (g) 


where  is  the  net  water  transport  from  anode  to  the 

cathode  through  the  electrolyte  membrane. 

In  the  cathode  catalyst  layer,  the  species  present  are  O2, 
N2,  and  H2O.  Nitrogen  does  not  react  in  the  cathode  cata¬ 
lyst  layer  and  the  amount  of  oxygen  consumed  and  water 
produced  can  be  calculated  with  Faraday’s  law: 


The  fuel  and  oxidant  undergo  electrochemical  reactions 
during  which  reactants  are  consumed  and  water  is  produced. 
The  amount  of  consumption  in  the  anode  and  cathode  catalyst 
layers  differ,  but  in  general  they  can  be  written  as: 

A  Nr=  J2  AA?r  (4) 

/=species 


where  the  summation  is  over  all  the  species  taking  part  in 
the  reaction.  For  the  anode  catalyst  layer,  the  species  present 
are  considered  to  be  H2  and  H2O.  H2  is  the  only  species 
consumed  in  the  anode  catalyst  layer.  The  rate  of  hydrogen 
molar  consumption  can  be  calculated  using  Faraday’s  law: 

AVrH2=-^^!i  (5) 

r  2  F 


where  /  is  the  current  density,  Ace n  is  the  active  area  of  a 
PEM  fuel  cell,  and  F  is  the  Faraday  constant. 


AN?2  =  I  (7) 

r  4  F 

ANf°  =  (8) 

r  2  F 

Water  molecules  dragged  from  the  anode  side  and  those 
produced  in  the  fuel  cell  are  assumed  to  enter  the  cathode 
gas  flow  channels.  Therefore,  the  molar  accumulation  in  the 
cathode  catalyst  layer  becomes: 


cathode  =  Aiy02  +  A  A^0  +  A^g  (9) 

As  the  water  content  of  the  cathode  side  flow  channels 
increases,  the  possibility  of  water  molecules  being  liquefied 
there  also  increases.  The  amount  of  liquid  water  present  is 
expressed  in  terms  of  the  thermodynamic  property  called 
quality,  and  the  thermodynamic  quality  in  each  segment 
of  the  cathode  side  is  determined  based  on  the  reactant 
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compositions,  local  temperature  and  pressure,  and  the  flow 
properties  are  modified  accordingly. 


2.2.2.  Conservation  of  energy 

In  order  to  satisfy  the  conservation  of  energy,  the  sum  of 
pressure  changes  around  each  of  the  loops,  i,  should  be  zero: 

^branch 

^  OijAPij  =  0  C*  =  1,2,3 . Moop)  (10) 

j=  i 

where  Ojj  is  a  sign  convention  representing  the  direction  of 
flow  in  the  branch  j  of  loop  i.  6j  j  is  considered  to  be  +1  when 
fluid  flows  in  clockwise  direction  and  —  1  if  the  direction  is 
reversed.  Total  pressure  drop  in  branch  j  of  loop  i  is  composed 
of  frictional,  APfjj,  and  local  or  minor  losses,  A Pmjj,  and 
can  be  written  in  the  form  of 


~  APf  Xj  +  AP] 


m  ,i,j 


(ID 


The  frictional  pressure  loss  can  be  expressed  by  Darcy- 
Weisbach  equation 


(12) 


where  Lfj  and  DhXj  are  the  branch  (or  control  volume) 
length  and  hydraulic  diameter,  pi  j  is  the  fluid  average  den¬ 
sity,  and  Vij  is  the  average  velocity  in  the  branch  j  of  the  loop 
i. 

The  friction  coefficient,  Cfjj,  is  a  function  of  the 
Reynolds  number  [5]: 


16  Ref  Ret  ,•  <  2000 

l,J  OJ  — 

0.079  Ref 14  Ret  ,  >  4000 

l,J  OJ  — 


(13) 


where  the  Reynolds  number  is  defined  based  on  the  hydraulic 
diameter: 


Pi,jVijDh,iJ 

JliJ 


(14) 


The  density  and  viscosity  of  the  fluid  are  calculated  incorpo¬ 
rating  the  flow  composition  and  operating  conditions  using 
the  equations  described  in  Baschuk  and  Li  [3].  A  linear  rela¬ 
tionship  is  used  to  calculate  friction  coefficient  in  the  transi¬ 
tion  regime  where  the  Reynolds  number  falls  between  2000 
and  4000. 

Minor  losses  arise  as  a  result  of  sudden  or  gradual  changes 
in  flow  directions  leading  to  the  formation  of  a  secondary 
recirculatory  fluid  motion.  Unlike  frictional  losses,  minor 
losses  are  local  and  occur  in  regions  where  the  recircula¬ 
tions  take  place.  The  relative  magnitude  of  minor  to  frictional 
losses  depends  primarily  on  the  hydraulic  network  geometry 
and  flow  conditions.  For  hydraulic  systems  with  long  ducts, 
minor  losses  may  be  neglected.  In  a  PEM  fuel  cell  stack,  on 
the  other  hand,  the  top  and  bottom  manifold  branches  are 
short,  the  change  in  the  flow  direction  is  sudden,  and  the 
number  of  turns  in  the  flow  channels  is  large;  as  a  result,  mi¬ 
nor  losses  may  become  major  and  should  be  included  in  the 


Table  1 

Branch,  confluence,  and  bend  head  loss  correlations  [6] 


-*¥r2 


ff 


p\V} 

APm,i-j  =  k\-j  ^  O'  =  2,  3),  k\—2 

-  10-“-1’-,‘(w)  +  i‘-2(w)! 

103  (V2\2  7000  fV3\2 

x  rT2  )  ’  1-3  “  IvJ 
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_  2100  /  Vi  \  2  _  7300  (V3\2 

-~R^\V2)  ,h~3  ~  ~R^  \V2) 

pV 2  7000 

APbend  —  2&i_3— —  ,0-3  =  — — 


analysis.  More  specifically,  these  so-called  minor  losses  can 
be  classified  into  three  groups:  (a)  minor  losses  due  to  the 
change  in  the  flow  direction  in  the  sharp  turns/bends  in  the 
gas  flow  channels,  and  the  minor  losses  associated  with  the 
(b)  confluence  and  (c)  branching  when  the  flow  channels  are 
either  charged  from  the  top  manifold  or  discharged  into  the 
bottom  manifold. 

Minor  losses  are  often  measured  experimentally  and 
expressed  by  empirical  correlations.  A  comprehensive  set 
of  formulations  for  laminar-to-laminar  and  turbulent-to- 
turbulent  flow  regimes  can  be  found  in  Jamison  and  Ville- 
monte  [6],  Blevins  [7]  and  Idelchik  [8].  For  fuel  cell  stack  ap¬ 
plications,  although  a  laminar  flow  prevails  in  the  flow  chan¬ 
nels  over  a  wide  range  of  operating  conditions,  a  combination 
of  turbulent-laminar  flow  regime  change  might  occur  at  the 
manifold-flow  channel  intersections  near  the  stack  inlet(s)  at 
large  current  densities.  The  minor  loss  coefficients  associated 
with  such  flow  regime  changes  have  not  yet  been  reported. 
Hence,  the  formulations  given  by  Jamison  and  Villemonte  [6] 
for  laminar  flow  as  listed  in  Table  1  are  used  in  the  present 
work.  The  locations  of  the  minor  losses  are  specified  by  “x” 
in  the  table.  It  is  seen  that  both  frictional  and  minor  losses 
depend  on  the  flow  distribution  in  the  stack. 


3.  Numerical  procedure 

Due  to  the  non-linear  relationship  between  the  pressure 
losses  and  mass  flow  rates  in  the  stack  manifolds  and  gas  flow 
channels,  the  mass  flow  rate  and  pressure  distributions  must 
be  calculated  using  an  iterative  technique.  The  flow  network 
solution  algorithm  devised  in  this  work  is  based  on  the  Hardy 
Cross  method  [9] ;  however,  it  has  been  extensively  modified 
to  account  for  the  variations  of  pressure  and  composition  and 
to  include  the  effect  of  minor  losses.  The  modified  algorithm 
requires  that  all  the  conservation  equations  be  written  in  terms 
of  total  molar  flow  rates.  This  approach  simplifies  tracking 
of  all  components  in  the  different  parts  of  the  stack. 

The  numerical  solution  begins  with  assigning  a  temporary 
flow  direction  to  the  flow  network.  The  sum  of  the  molar  flow 
rates  from  all  the  inlets  (with  known  compositions)  is  then 
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N  inlet.i 


Fig.  4.  Schematic  illustration  for  nodal  mass  conservation. 


divided  uniformly  among  the  upstream  sections  of  all  the  gas 
flow  channels.  The  molar  flow  rates  and  compositions  in  the 
downstream  sections  of  the  flow  channels  are  calculated  by 
subtracting  the  consumed  reactants  and  adding  the  produced 
components.  Water  transport  through  the  membrane  due  to 
the  combined  effects  of  electroosmotic  drag  and  back  diffu¬ 
sion  is  also  considered  (Eqs.  (6)  and  (9)).  The  next  step  is  to 
use  the  conservation  of  mole  equation  and  calculate  the  molar 
flow  rates  in  the  top  and  bottom  manifolds  based  on  the  prede¬ 
fined  directions  and  the  assumed  molar  flow  rates  in  the  flow 
channels.  Fig.  4  illustrate  three  nodes  in  the  top,  middle  (reac¬ 
tion  site),  and  bottom  of  a  typical  flow  channel.  Application 
of  the  conservation  of  mole  equation  to  these  nodes  result  in 
total  molar  flow  rates,  Nt,  in  branches  1  and  4  of  the  ith  loop 
and  branch  3  of  the  (/  —  l)th  loop,  respectively,  as  follows: 

vr  @i— \,\Ni— 11  +  Wniet ,i  ~  ft— 1,2^7— 1,2  /1 

Ni,i  =  - 1 - : - — - - : - -  (15) 

“i,  1 


Ni- 1,3  = 


0: 


i- 1,3 


(16) 


Ni,  4  = 


Oj-l,  4^7-1, 4  ^outlet,/ 

OiA 


(17) 


where  Mniet,/  and  7V0Utiet,f  are  the  total  molar  flow  rates 
entering  and  leaving  the  stack  at  that  particular  location. 
A  negative  molar  flow  rate  (in  any  of  the  above  equations) 
indicates  that  the  assumed  flow  direction  in  that  branch  is 
incorrect  and  must  be  reversed. 

The  conservation  of  energy  equation  must  also  be  satisfied 
by  the  assumed  molar  flow  rates.  Eq.  (10)  can  be  expressed 
in  terms  of  molar  flow  rates  as  follows: 


^branch 

*Pi  =  E  MauC  +  +  '-3./,/>iv"yui  =  o 

j= 1 

(i=  1,2,3,...,  Moop)  (18) 


where  the  frictional  and  minor  losses  have  been  expressed  in 
the  form  r Nn .  For  every  branch  r\  is  the  coefficient  associ¬ 
ated  with  the  frictional  loss,  r 2  is  the  coefficient  associated 
with  the  branch  or  confluence  loss,  and  r 3  represents  the  co¬ 
efficient  for  the  losses  in  the  bends  of  gas  flow  channels.  The 
r  and  n  parameters  are  obtained  by  rewriting  Eq.  (12)  and  the 
correlations  listed  in  Table  1  on  the  molar  flow  rate  basis.  For 
instance,  r\jj  can  be  expressed  as: 


f  0.0407 LijfajRT 


ruj—  < 


n 4  p.  . 


(laminar) 


1.3560  x  10"6L;  upf  p?-]5(RT)L15 

n4.75  pi. 75 


(turbulent) 


(19) 


where  R  is  the  gas  constant.  T  is  the  stack  operating  tempera¬ 
ture  and  considered  to  be  constant  throughout  the  stack.  Pi  j 
represents  the  absolute  pressure  of  the  branch.  The  value  of 
n  is  equal  to  1  for  laminar  flow  and  1.75  for  turbulent  flow.  It 
is  seen  that  the  value  of  r\  depends  exclusively  on  the  con¬ 
ditions  of  the  branch  it  is  calculated  for  and  independent  of 
other  branches. 

Minor  losses  can  be  expressed  similarly;  however,  un¬ 
like  frictional  losses,  the  minor  loss  of  a  branch  depends 
not  only  on  the  flow  conditions  of  the  branch  but  also  is 
influenced  by  the  conditions  of  the  neighboring  branches. 
Therefore,  extra  care  must  be  taken  to  include  the  influence 
of  neighboring  branches  in  the  calculations.  Table  2  is  equiv¬ 
alent  to  Table  1  in  which  all  possible  minor  losses  are  pre¬ 
sented  on  molar  flow  rate  basis.  As  indicated,  the  power  of 


Table  2 

Branch,  confluence,  and  head  loss  correlations  rewritten  in  terms  of  the  molar  flow  rate,  otherwise  identical  with  those  shown  in  Table  1 
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the  molar  flow  rates  for  all  different  minor  losses  is  equal  to  1 

(*2  =  I)- 

Although  the  initially  assumed/calculated  molar  flow  rates 
satisfy  the  conservation  of  mass,  the  conservation  of  energy 
may  not  be  satisfied.  Hence,  the  corrected  molar  flow  rates, 
Njj,  can  be  related  to  the  assumed  flow  rates,  Nijto,  as: 

Nij  =  NiJt  o  +  A  Ni  (20) 

where  A  A;  is  the  correction. 

For  each  branch,  total  pressure  loss  can  be  expressed  as: 


A  Pij  = 


rujWj,  o  +  A  V  ' 

+(r2  JJ  +  r3,i,j)(.UijQ  +  A  Ni) 

/  V  1  i  i  V  1  i  /  1  a  V  t 

rU,/Ni,j,0  +nU,jNi,j,6  AiVii- 

+(r2,(,;  +  r3jj)(Njj:o  +  ANj) 


if  A  Ni  is  small  compared  with  all  terms  of  the  series 

after  the  second  can  be  dropped.  The  conservation  of  energy 
principle  requires  that  for  every  loop  i  in  the  flow  network 


^branch 

E  'V, A 

7=1 

^branch 

=  ^  @i,j(r\JJ Nif j' Q  +  j, 0 ) 

7=1 


+  AA; 


^branch 


^  j j’ Q7  +  (r2,i,j  +  r3,i,j)) 


7=1 


—  0  (/  —  1,  2,  3, ... ,  A^ioop) 


(21) 


in  which  ANj  has  been  taken  out  of  the  summation  because 
it  is  the  same  for  all  branches  j  =  1,  2,  3  ... ,  A^ranch  in  the 
loop  i.  Eq.  (21)  is  solved  for  ANj  in  each  loop  in  the  network 
as 


Ejr  o'  +  +  rsjjWij,  o) 


AM-  =  - 


S  /=jinth  ®i,  j(n  u U  j r  1 , i,  j fy,  j, o  +  (r2, i,  j  +  r3 , i,  j )) 

(7  —  1,  2,  3,  ... ,  A^i00p)  (22) 


AT^l.U-l 


With  the  initial  guesses  for  the  flow  distributions,  it  is  ex¬ 
pected  to  have  nonzero  molar  flow  rate  corrections.  Hence, 
the  numerical  solution  proceeds  as  follows:  the  values  of 

AiV’s  are  calculated  for  all  the  loops  in  the  flow  network 
and  the  molar  flow  rates  are  modified  accordingly.  For  the 
branches  shared  by  two  successive  loops,  the  correcting  mo¬ 
lar  flow  rates  from  both  loops  must  be  accounted  for.  The 
iterations  are  continued  until  convergence  is  reached.  For  the 
present  study  convergence  is  considered  to  have  been  reached 
if  the  maximum  of  A  P’s  for  all  the  loops  is  less  than  1  Pa,  as 
shown  below: 


Maximum(|  AP/|,  /  =  1,2,3,...,  Moops)  E  1  Pa  (23) 


The  solution  of  the  stack  flow  network  model  provides  the 
molar  flow  rate  distributions,  composition  and  pressure  in  the 
flow  channels.  These  values  are  used  by  the  single  cell  model 
to  calculate  Pceii  (and  ^cp)  for  each  cell.  Finally,  the  stack 
voltage  is  obtained  by  Eq.  (1). 


4.  Results  and  discussion 

The  input  parameters  for  the  fuel  cell  and  stack  flow  model 
are  classified  as  operating  and  design  parameters.  The  design 
parameters  are  the  fuel  cell  size,  stack  manifold  and  flow 
channel  dimensions  and  configuration.  Operating  parame¬ 
ters  include  the  stack  current  density,  temperature,  pressure, 
stoichiometry  and  the  reactant  composition  at  the  stack  inlet. 
Table  3  lists  a  summary  of  the  operating  and  relevant  design 
parameters  for  various  cell  and  stack  components  such  as  the 
bipolar  plate,  electrode  backing,  catalyst  layer  and  polymer 
electrolyte  membrane. 

Fig.  5  shows  the  effect  of  minor  losses  on  the  pressure 
distributions  along  the  top  and  bottom  manifolds  of  a  U- 
configuration  PEM  fuel  cell  stack  as  shown  in  Table  5  (Run 
#1  without  minor  losses  and  Run  #2  with  minor  losses).  The 
stack  consists  of  51  cells  and  operates  at  5000  Am-2.  Fig. 
5(a)  and  (b)  shows  the  pressure  distributions  for  the  anode 
and  cathode  sides  of  the  stack  without  considering  the  mi¬ 
nor  losses.  The  inlet  pressures  are  at  the  highest  and  drop 

Table  3 


Parameters  and  properties  used  in  the  present  PEM  fuel  cell  stack  simulations 


Component 

Parameter 

Value 

Bipolar  plate 

W 

0.12m 

H 

0.12m 

l 

0.3-4. 8  m 

d\\ 

5  x  10-4  to  2  x  10-3  m 

Nc 

1-10 

Nt 

3-40 

Electrode  backing3 

nbulk 

A'e 

6  x  10~5  Qm 

<$e 

2.5  x  10~4  m 

0e 

0.4 

Catalyst  layer3 

2.0465  x  10“5  m 

rapt 

0.004  kg  m2 

/pt 

0.2 

U 

0.9 

Ks 

727008  m-1 

Polymer  electrolyte  membrane3 

5m 

1.64  x  10-4  m 

7.18  x  10-20  m2 

KP 

1.8  x  10-18  m2 

CH+ 

1200  mol  m-3 

£>h+ 

4.5  x  10~9  m2  s~2 

Stack 

U cells 

51 

U inlets 

1-3 

T 

353  K 

Pout 

1  atm 

J 

1000-10000  Am"2 

<Sa 

1.2 

<SC 

2.0 

a  Parameters  associated  with  the  single  cell  model  [3]. 
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Fig.  5.  Effect  of  minor  losses  on  the  pressure  distributions  in  a  U-configuration  stack:  (a  and  b)  anode  and  cathode  sides’  pressure  distributions  without  the 
minor  losses,  and  (c  and  d)  anode  and  cathode  sides’  pressure  distributions  with  the  minor  losses  (Runs  #1  and  #2:  D h  =  0.01  m,  <4  =  0.002  m,  /  =  0.3  m, 
Nc  —  10,  Nt  =  3  and  J  =  5000  A m~2). 


gradually  along  the  manifolds  as  flows  are  discharged  into 
the  flow  channels.  Significant  pressure  drops  occur  within 
the  flow  channels  as  reactants  pass  through  the  tiny  channels. 
The  unreacted  components  are  then  exited  into  the  bottom 
manifold,  collected  channel  by  channel,  and  exit  the  stack 
at  atmospheric  pressure.  Fig.  5(c)  and  (d)  shows  the  effects 
of  minor  losses  on  the  pressure  distributions  in  the  anode 
and  cathode  sides.  With  the  minor  losses  present,  larger  inlet 
pressures  are  required  to  drive  flows  in  the  stack.  Numerical 
results  indicate  that  the  minor-to-frictional  loss  ratios  are  in 
the  order  of  4  for  both  anode  and  cathode  sides  of  the  stack. 
Since  the  flow  and  pressure  distributions  within  the  stack  are 
strongly  coupled,  it  is  therefore  expected  that  minor  losses 
play  an  important  role  in  the  stack  performance. 

Essentially,  a  large  pressure  drop  in  the  manifolds  may 
lead  to  uneven  gas  flow  distribution  in  the  stack.  Uneven  flow 
distribution  causes  a  low  reactant  concentration  in  some  of 
the  cells  and  degrades  the  performance  of  these  cells,  hence 
the  overall  performance  of  the  stack.  On  one  hand,  if  the 
pressure  drop  in  the  flow  channels  on  the  bipolar  plates  is 
large,  the  manifold  behaves  like  a  plenum  leading  to  a  more 


uniform  reactant  and  hence  stack  voltage  distributions.  The 
large  pressure  drop,  on  the  other  hand,  puts  extra  load  on 
the  air  compressor.  Therefore,  the  ratio  of  the  pressure  loss 
in  the  flow  channels  to  that  of  the  manifold,  herein  after  de¬ 
noted  by  a,  is  a  good  indication  of  how  the  flow  is  distributed 
throughout  the  stack  and  how  the  stack  performs  in  general. 
The  value  of  a  is  particularly  important  at  the  oxidant/cathode 
side  since  the  resistance  of  flow  of  hydrogen  is  very  low  due 
to  the  low  flow  rate  involved  and  therefore  hydrogen  pressure 
losses  are  usually  insignificant. 

Table  4  summarizes  the  pressure  data  from  Fig.  5  and 
lists  the  a  values  for  comparison.  As  indicated  in  this  table, 
without  considering  minor  losses,  this  ratio  is  much  larger 
for  the  anode  side  than  that  of  the  cathode  side.  Hence  a 
more  uniform  flow  distribution  is  expected  to  prevail  in  the 
anode  side  of  the  stack.  However,  when  the  minor  losses  are 
included,  the  ratio  decreases  considerably  for  the  anode  side 
and  only  a  little  for  the  cathode  side.  As  a  result,  the  presence 
of  minor  losses  degrades  the  flow  distribution  in  the  anode 
side  without  affecting  significantly  the  flow  distribution  in 
the  cathode  side.  This  is  illustrated  in  Fig.  6(a)  and  (b)  where 


Table  4 


Relative  magnitudes  of  the  frictional  and  minor  losses  in  different  sections  of  a  U-configuration  stack  operating  at  5000  Am 


A  P  (Pa) 

Manifolds 

Flow  channels 

Stack 

Frictional 

Minor 

Total 

Frictional 

Minor 

Total 

af 

at 

Anode 

30 

591 

621 

184 

846 

1030 

6.17 

1.6 

Cathode 

360 

1472 

1833 

939 

3735 

4674 

2.60 

2.5 
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Table  5 

Stack  performance  for  a  variety  of  stack  configurations  operating  at  5000  Am 


Run  no.  Stack  configuration  AP  (Pa) 


la 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


U-Configuration 


1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

II 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

II 

1 1 

1 1 

1 1 

1 1 

II 

1 1 

1 1 

ll-Configuration 


1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 
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Fig.  6.  Effect  of  minor  losses  on  the  stoichiometry  distributions  in  a  U-configuration  stack:  (a)  anode  stoichiometry  and  (b)  cathode  stoichiometry  (Runs  #1, 
#2:  Dh  =  0.01  m,  <4  =  0.002m,  /  =  0.3m,  Nc  =  10,  Nt  =  3  and  J  —  5000  Am-2). 
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cell-to-cell  stoichiometry  distributions  are  plotted  for  anode 
and  cathode  sides,  respectively. 

Fig.  7  displays  the  fuel  and  oxidant  molar  flow  rate  and 
concentration  distributions  in  the  top  and  bottom  manifolds 
of  the  anode  and  cathode  sides  of  the  stack.  The  effect  of 
minor  losses  on  the  distributions  are  also  included  in  the  fig¬ 
ure.  As  shown  in  Fig.  7(a)  and  (b),  the  H2  and  O2  molar  flow 
rates  are  maximum  at  the  stack  inlet  and  decrease  along  the 
top  manifold  as  the  flow  channels  are  fed  with  the  reactants. 
A  similar  trend  is  observed  in  the  bottom  manifold  as  unre¬ 
acted  components  along  with  the  product  water  are  collected 
from  the  flow  channels  towards  the  stack  outlet.  The  con¬ 
centration  profiles  calculated  based  on  the  control  volumes’ 
molar  volumetric  flow  rates  are  shown  for  H2  and  O2  in  Fig. 
7(c)  and  (d),  respectively.  It  is  evident  from  these  figures  that 
the  effect  of  minor  losses  is  more  significant  for  the  anode 
side  of  the  stack  as  just  pointed  out  earlier. 

The  effect  of  minor  losses  on  the  fuel  and  oxidant  molar 
flow  rates  and  concentration  distributions  in  the  upstream  and 
downstream  flow  channels  are  shown  in  Fig.  8.  Again,  it  is 
seen  that  the  presence  of  minor  losses  has  an  adverse  effect  on 
the  flow  distribution  in  the  anode  side  and  minimal  influence 
on  the  distributions  in  the  cathode  side. 

The  logarithmic  average  of  reactant  concentrations  along 
with  the  pressure  distributions  in  the  stack  flow  channels  are 
used  in  the  single  cell  model  to  calculate  stack  voltage  distri¬ 
butions  for  the  U-configuration.  Stack  temperature  is  consid¬ 
ered  to  be  constant  at  80  °C.  Fig.  9  shows  how  the  resulting 
voltage  varies  from  one  cell  to  another.  As  indicated  in  this 
plot,  the  cell  voltage  is  maximum  near  the  stack  inlet  and 
decreases  as  the  distance  from  the  inlet  increases.  The  pres¬ 
ence  of  the  minor  losses  in  the  system  tends  to  modify  the 
voltage  distribution  and  increase  the  cell-to-cell  voltage  vari¬ 
ation.  The  cell-to-cell  voltage  variation  can  be  quantified  by 
the  voltage  spread  of  the  stack  defined  as 
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Fig.  10.  Effect  of  inlet  location  on  the  (a)  anode  stoichiometry,  (b)  cathode 
stoichiometry,  and  (c)  cell-to-cell  voltage  distributions  (Runs  #2  to  #7:  D h  = 
0.01  m,  <4  —  0.002m,  /  —  0.3m,  Nc  =  10,  Nt  =  3  and  J  =  5000  Am-2). 


Fig.  9.  Cell-to-cell  voltage  distributions  in  the  U-configuration  stack  (Runs 
#1  and  #2:  D h  =  0.01  m,  d h  =  0.002  m,  /  =  0.3  m,  Nc  =  10,  Nt  =  3  and 
/  =  5000  Am-2). 


where  E™f  and  are  the  maximum  and  minimum  cell 
voltages,  respectively,  within  the  stack.  For  the  present  U- 
configuration  and  operating  conditions,  the  voltage  spread 
is  increased  by  about  10%  from  4.2%  without  considering 
system  minor  losses  to  4.6%  when  minor  losses  are  taken 
into  account.  This  figure  also  shows  that  a  slightly  larger 
stack  voltage  (about  0.3%)  is  achieved  when  minor  losses 
are  taken  into  account. 

High  voltage  spread  is  an  indication  of  stack  poor  perfor¬ 
mance.  In  practice,  a  low  voltage  spread,  uniform  cell-to-cell 
performance,  and  low  energy  loss  due  to  the  pressure  drops 
are  desirable  in  order  to  maximize  stack  performance.  There¬ 
fore,  it  is  very  helpful  to  examine  other  stack  configurations 
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cathode  side  of  a  symmetric  double  inlet  stack  (Run  #9:  D h  =  0.01  m,  d h  =  0.002  m,  /  =  0.3  m,  Nc  =  10,  Nt  =  3  and  J  =  5000  Am-2). 


to  study  the  effect  of  inlet  and  outlet  locations,  the  effect  of 
number  of  inlets,  and  bipolar  flow  channel  design  conditions 
on  the  PEM  fuel  cell  stack  performance. 

With  minor  losses  considered,  the  stack  inlet  is  moved 
from  the  present  location  (U-configuration)  towards  the  stack 
other  end-plate  (Z-configuration)  and  the  resulting  anode  and 
cathode  side  stoichiometries  and  cell-to-cell  voltage  distribu¬ 
tions  are  calculated.  Fig.  10(a)  and  (b)  depicts  variations  in 
the  anode  and  cathode  stoichiometries  as  functions  of  the  in¬ 
let  location.  The  resulting  voltage  distributions  are  displayed 
in  Fig.  10(c).  It  is  clear  from  these  figures  that  all  the  dis¬ 


tributions  are  improved  by  moving  the  inlet  away  from  the 
outlet.  When  inlet  is  placed  close  to,  or  at  the  opposite  side 
of  the  outlet  (e.g.  Z-configuration),  the  PEM  fuel  cell  stack 
operates  with  the  minimum  voltage  spread.  This  is  in  agree¬ 
ment  with  the  results  recently  reported  by  Baschuk  and  Fi 
[4]  in  the  absence  of  minor  losses.  Other  stack  configurations 
with  variant  inlet/outlet  topologies  were  examined  and  the 
numerical  results  are  summarized  in  Table  5.  In  this  table,  in 
addition  to  the  stack  voltage  data,  the  average  pressure  drops 
for  the  cathode  side  of  the  stack  and  representative  compres¬ 
sor  power  requirements  for  the  oxidant  supply,  Pc ,  are  also 


Fig.  12.  Pressure  distributions  along  the  top  and  bottom  manifolds  of  the  symmetric  double  inlet  stack:  (a)  anode  side,  and  (b)  cathode  side  (Run  #9:  D h  =  0.01  m, 
dh  =  0.002m,  /  —  0.3  m,  Nc  =  10,  Nt  =  3  and  J  —  5000  Am-2). 
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Fig.  13.  Variations  of  (a)  Vh2  and  (c)  Ch0  in  the  upstream  and  downstream  flow  channels  of  the  anode  side  and  (b)  Non  and  (d)  Co2  in  the  upstream  and 
downstream  flow  channels  of  the  cathode  side  of  a  symmetric  double  inlet  stack  (Run  #9:  D h  =  0.01  m,  d h  =  0.002  m,  /  =  0.3  m,  Nc  =  10,  Nt  =  3  and 
/  =  5000  Am-2). 


calculated  for  different  configurations.  It  is  evident  from  this 
table  that  a  two-inlet  configuration  (Run  #9)  leads  to  the  low¬ 
est  voltage  spread.  Hence  this  configuration  has  been  chosen 
for  further  investigations. 

Fig.  11(a)  and  (b)  depicts  variations  of  fuel  and  oxidant 
molar  flow  rates  in  the  manifolds  of  the  anode  and  cathode 
sides,  respectively.  The  reactants  molar  flow  rates  are  de¬ 
creased  linearly  from  the  inlets  (located  in  both  end-plates) 
towards  the  manifold  midpoint  indicating  even  molar  distri¬ 
butions  in  the  flow  channels.  Electrochemical  reaction  occurs 
at  the  reaction  site  of  each  cell,  and  the  unreacted  components 
and  the  reaction  product  are  collected  from  the  flow  channels 
and  accumulated  as  they  move  towards  the  stack  outlet(s)  at 
the  center  of  the  stack.  The  corresponding  H2  and  O2  con¬ 
centration  profiles  calculated  in  the  control  volumes  along 
the  manifolds  are  displayed  in  Fig.  11(c)  and  (d). 

The  pressure  variation  in  the  stack  plays  an  important  role 
in  the  flow  distributions  in  the  flow  channels.  Fig.  12(a)  and 
(b)  illustrate  how  the  pressure  varies  along  the  anode  and  cath¬ 
ode  top  and  bottom  manifolds.  Careful  examination  of  pres¬ 
sure  data  indicates  that  the  flow  channel-to-manifold  pressure 
drop  ratios,  a ,  is  about  6.2  and  9.2  for  the  anode  and  cath¬ 
ode  sides,  respectively.  The  large  a  values  for  the  symmetric 
double-inlet  stack  configuration  (as  compared  with  1.6  and 
2.5  for  U-configuration  in  Table  4)  is  a  good  indication  of  uni¬ 
formity  in  the  flow  distribution  and  verifies  the  linear  drop 


in  the  molar  flow  rates  in  the  anode  and  cathode  manifolds 
as  shown  in  Fig.  11(a)  and  (b).  The  corresponding  fuel  and 
oxidant  molar  flow  rates  and  concentration  distributions  are 
shown  in  Fig.  13.  As  seen  in  these  figures,  the  symmetric  dou¬ 
ble  inlet  configuration  forces  the  fuel  and  oxidant  molar  flow 
rates  and  concentrations  to  increase  in  the  stack  midpoints 
as  shown  in  Fig.  13(a)  and  (b).  As  a  result,  more  uniform 
fuel  and  oxidant  concentrations  are  formed  in  the  stack  as 
indicated  in  Fig.  13(c)  and  (d),  respectively. 

The  sizing  optimization  of  the  bipolar  flow  channels  for 
a  given  basic  geometry  is  a  complex  task  in  which  sev¬ 
eral  requirements  must  be  balanced.  The  issue  generally 
requires  more  consideration  on  the  cathode  side.  The  op¬ 
timization  of  channels  is  not  very  thoroughly  discussed  in 
the  literature.  Some  general  issues,  however,  can  be  out¬ 
lined.  In  the  present  study,  the  effects  of  the  manifold  and 
flow  channel  design  conditions  are  studied  for  the  two-inlet 
stack  and  bipolar  plate  configuration  illustrated  previously  in 
Fig.  2. 

Fig.  14(a)  shows  the  strong  effect  of  manifolds’  hydraulic 
diameter  on  the  cathode  pressure  loss  ratio.  For  D h  >  0.01  m 
the  ratio  (a)  increases  considerably  such  that  the  top  mani¬ 
fold  behaves  like  a  plenum,  and  distributes  reactants  uni¬ 
formly  into  the  flow  channels.  As  a  result,  cell-to-cell  volt¬ 
age  distribution  becomes  more  uniform  as  displayed  in  Fig. 
14(b).  However,  if  D h  <  0.005  m,  the  pressure  drop  in  the 
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Fig.  14.  Variation  of  (a)  a,  and  (b)  voltage  spread  as  functions  of  manifolds 
hydraulic  diameters  in  a  symmetric  double  inlet  stack  (dh  =  0.002  m,  /  = 
0.3m,  Nc  =  10,  Nt  =  3  and  J  =  5000  Am-2). 
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top  manifold  is  comparable  with  (or  even  less  than)  that  of 
flow  channels,  leading  to  a  larger  voltage  spread.  Numerical 
results  indicated  that  enlarging  the  manifold  diameter  does 
not  change  significantly  the  energy  requirement  of  the  stack 
(e.g.  compressor  power)  hence  it  can  be  concluded  that  a 
hydraulic  diameter  between  0.01  m  and  0.02  m  is  optimum 
for  which  the  least  voltage  spread  is  achieved  at  minimum 
compressor  power  requirement. 

Numerical  results  have  also  provided  some  guidelines  for 
bipolar  plate  design.  The  effect  of  number  of  flow  channels 
per  bipolar  plate  of  the  cathode  side  on  the  stack  perfor- 


Table  6 


Effect  of  bipolar  design  parameters  on  the  performance  of  a  double-inlet 
stack  operating  at  5000  Am-2 


Manifolds, 
Dh  (m) 

Bipolar  Plate 

Stack 

Nc 

Nt 

/(m) 

dh  (m) 

a 

Se  (%) 

Estack 

(V) 

Pc  (W) 
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Fig.  15.  Variations  of  (a)  voltage  spread,  (b)  stack  voltage  and  (c)  estimated 
power  requirement  for  oxidant  supply  as  functions  of  current  density  for  a 
symmetric  double  inlet  stack  (Dh  =  0.01  m,  d h  =  0.002  m,  /  =  0.3  m,  and 
Nc  =  10,  At  =  3). 

mance  are  indicated  in  Table  6.  The  flow  channel  diameter, 
dh,  length,  /c,  number  of  flow  channels,  Nc  and  number  of 
turns,  Nt  are  all  related  (Fig.  2).  If  the  number  of  flow  chan¬ 
nels  is  decreased,  the  flow  channel  length(s)  and  the  number 
of  turns  must  be  increased  (to  cover  the  same  transport  area) 
both  of  which  augment  the  minor  and  frictional  losses  in  the 
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bipolar  plate.  Although  the  increased  pressure  loss  in  the  gas 
flow  channel  improves  the  stack  voltage  distribution  and  re¬ 
duces  voltage  spread  and  helps  to  remove  water  droplets  in 
the  stack,  its  obvious  drawback  is  that  the  system  efficiency 
is  reduced  because  the  power  consumption  of  the  air  supply 
is  increased  as  indicated  in  the  table.  Numerical  results  indi¬ 
cate  that  for  small  number  of  flow  channels  the  stack  voltage 
spread  remains  close  to  zero. 

The  stack  performance  is  also  affected  by  the  cell’s  cur¬ 
rent  density.  Fig.  15(a)  and  (b)  shows  that  the  voltage  spread 
increases  and  the  stack  voltage  decreases  as  the  current  den¬ 
sity  is  increased.  The  system  pressure  drop  and  the  power 
requirement  increase  at  higher  current  densities,  since  larger 
reactant  molar  flow  rates  are  now  fed  into  the  stack.  Variation 
of  the  power  requirement  versus  the  current  density  is  shown 
in  Fig.  15(c). 


5.  Conclusions 

The  performance  of  a  PEM  fuel  cell  stack  consisting  of 
51  active  cells,  with  various  inlet-outlet  configurations,  was 
studied  under  various  design  and  operating  conditions.  Mo¬ 
lar  flow  rate  and  pressure  distributions  in  the  stack  manifolds 
and  gas  flow  channels  are  determined  using  a  modified  ver¬ 
sion  of  Hardy  Cross  method.  Minor  losses  due  to  the  sharp 
turns  in  the  flow  channels  and  confluence/branching  flows  at 
the  manifold-flow  channel  intersections  are  considered.  The 
results  of  the  flow  network  analysis  were  used  in  a  single 
cell  model  previously  developed  to  evaluate  the  stack  perfor¬ 
mance.  From  the  numerical  results,  the  following  conclusions 
can  be  made.  Minor  losses  increase  the  stack  operating  pres¬ 
sure  and  the  compression  load  and  alter  the  voltage  spread  in 
the  stack.  The  inlet-outlet  configuration  plays  an  important 


role  in  the  stack  performance.  Although  a  Z-configuration 
stack  is  preferable  to  U-configuration,  simulation  results  have 
revealed  that  a  symmetric  double  inlet-single  outlet  topology 
provides  excellent  performance  with  reasonably  low  com¬ 
pressor  requirement  demand  and  minimum  voltage  spread. 
The  stack  performance  can  also  be  quantified  by  the  bipo¬ 
lar  plate-to-manifold  pressure  loss  ratio,  a.  For  a  values  of 
greater  than  10,  the  stack  voltage  spread  of  less  than  1  %  can  be 
achieved.  The  size  and  the  number  of  flow  channels  carved  on 
the  bipolar  plate  are  also  very  important.  Flow  channels  with 
smaller  cross-sectional  area  and  longer  lengths  increase  the 
compressor  power  demands  considerably  without  improving 
the  stack  voltage  spread  appreciably. 
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